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We consider the dynamics of an isolated quantum many-body system after a sudden change of one 
control parameter, focusing on the statistics of the work done. We establish a connection between the 
generating function of the distribution of the work and the partition function of a classical system in 
a film geometry. Using this connection, we first show that the scaling of the fidelity susceptibilities 
close to a quantum phase transition can be understood in terms of the critical behavior of the 
excess entropy and specific heat in the classical model. Remarkably, we show that the statistics of 
the work close to the threshold and to criticality is connected to the so-called critical Casimir free 
energy which is responsible for the interaction between the boundaries of the classical system. On 
the basis of this relation, we highlight the emerging universal features of the statistics of the work. 
Our findings are exemplified for a global quench of the transverse Ising chain at zero temperature. 

PACS numbers: 05.70.Ln, 03.65.Yz. 



Universality, i.e., the insensitivity to microscopic de- 
tails, is one of the key discoveries of modern condensed 
matter and statistical physics. We definitely owe to this 
concept the possibility to describe the low-temperature 
behavior of strongly correlated quantum systems in terms 
of simple effective models which account for just a few rel- 
evant interactions and symmetries. This is of crucial im- 
portance in many contexts, among others in the study of 
the low-temperature physics of real materials, where the 
overwhelming number of details to be accounted for in an 
ab initio description would definitely hinder all progress. 

Notwithstanding the insight into universality at equi- 
librium, the investigation of its emergence in the dynam- 
ics of quantum many-body systems started only recently, 
inspired by a series of ground breaking experiments in the 
context of cold atoms in optical lattices 0] ■ These sys- 
tems constitute a sort of ideal nonequilibrium laboratory: 
long coherence times and an exquisite control on the sys- 
tem made it possible to observe interesting phenomena, 
such as the collapse and revival cycles of a system taken 
abruptly across the Mott-superfiuid phase boundary [l| . 
A common protocol employed to take the system out 
of equilibrium is the so-called quantum quench, i.e., a 
rapid or slow change in time of one of the parameters 
controlling the system Hamiltonian A great deal of 
recent research focusing on these protocols has uncovered 
various examples of universal features, from the scaling 
of the energy density for slow quenches across quantum 
critical points Q to the the asymptotics of correlators af- 
ter a quench at a conformally invariant quantum critical 
point 0- 

In the light of thermodynamics, quantum quenches are 
essentially thermodynamic transformations 0, Q , which 
are naturally dscribed by three standard thermodynamic 
quantities: change in entropy, transferred heat [5|, and 
work W done on the system [6|. Focusing on the latter, it 



is well know both in classical and quantum physics that 
W is in general a random variable, fiuctuating in differ- 
ent realizations of the same protocol and described by 
a probability distribution P(W) 0, Q. The purpose of 
the present study is to investigate the emergence of uni- 
versal features in the statistics P{W) of the work done 
on the system by changing abruptly one of its parame- 
ters @, Q- In order to extract all universal aspects of 
P(W) we connect the generating function of the work 
distribution to the partition function of a classical sys- 
tem in a film geometry. We first use this mapping to 
gain insight into the scaling of the fidelity susceptibili- 
ties close to quantum critical points [loj . Afterwards, we 
discuss the universality of the statistics of the work W 
close to the threshold by uncovering its connection to the 
universal critical Casimir force (see, e.g., Ref. [ll[) which 



describes the interaction between the boundaries in the 
classical model. In order to substantiate our general find- 
ings we focus on the example of a quench of the transverse 
field in the quantum Ising chain close to criticality. 

Virtually all instances of quantum quenches consid- 
ered so far in the literature consist of the change of a 
parameter 5 of a many-body Hamiltonian H{g) from 
an initial value gi to a final one gf. In this set- 
ting, the work W done on this closed system can be 
quantified measuring the energy before and after the 
quench. If the quench starts in the ground state \'^o{gi)) 
of the initial Hamiltonian, the probability distribution 
of the work W is P{W) = En>o HW-[Enigf) - 
Eoig,)])\{^ 0(91)1-^ nig f))\^ where |*„(ff/)) are the eigen- 
states of energy En{gf) of the final Hamiltonian and P 
vanishes for W < AEq = Eoigj) — E^igi). The char- 
acteristic function G{t) 
takes the form 



^_^AWe''^*P{W) of P{W) 



Git) = {^oig^W'^''^'e-''"(Sf)t\^oig^)), (1) 
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where here t is the variable conjugate toW |G(i)P 
is known as Loschmidt echo and has been studied in var- 
ious contexts ranging from quantum chaos and quantum 
information to strongly correlated systems [l^ . 

Essentially, the characteristic function G{t) can be seen 
as a partition function after a proper analytic continua- 
tion to real "times" t. This simple, yet crucial observa- 
tion yields insight into the universal features of P{W). 
In order to see this, consider Eq. and perform an 
analytic continuation to imaginary "time" t — > —iR. We 
obtain 



G{R) = e 



_ ^-AEoR 



X Z{R), where 



(2) 



For a d-dimcnsional quantum system possessing a + 1- 
dimensional classical correspondent, Z{R) is the parti- 
tion function of the latter on a "strip" of thickness R 
with two boundaries described by the boundary states 
|^o(5i))- If the system has linear extension L it is now 
possible to separate three contributions to the free energy 
-lnG(i?) = L'i[i?x/b+2/,+/,(i?)], where /fa = ^E^/L'^ 
is the bulk free energy density and fs is the surface free 
energy associated to each of the two (identical) bound- 
aries of the film. The remaining contribution fc{R) repre- 
sents an effective interaction between these two bound- 
aries, which generically decays to zero at large separa- 
tions R. In particular, one can easily identify 



Z{R) 



(3) 



Upon approaching a critical point, fc{R) becomes univer- 
sal and causes the so-called critical Casimir effect [llj. 

Relating quantum quenches to the statistical physics 
of classical confined systems makes it possible to gain 
insight into the scaling of the fidelity susceptibility, a 
quantity recently introduced for the characterization of 
quantum phase transitions [lo| . The fidelity J-(gi,gf) = 
|('J'o(5i)l^o(5/))|i determines, inter alia, the weight of 
the peak of P{W) at the threshold W = AEf^ In- 
deed, if the quantum system has a critical point at 
g = gc, the generalized susceptibilities Xn[gi,gf) = 
—L~'^dg^ \TiIF{gi,gf) develop non-analytic behaviors 
and in particular the so-called fidelity susceptibility xi 
scales as 



X2(5,5) - \9~ 9c\' 



(4) 



This scaling can be straightforwardly understood in 
terms of the statistical mechanics of the correspond- 
ing d -\- 1-dimensional classical systems with boundaries, 
for which g acts as a temperature-like variable control- 
ling the distance from the bulk critical point. In fact, 
by inserting in Eq. ^ a resolution of the identity in 
terms of the eigenstates of H{gf) and using Eq. ^ for 
R — >■ -|-oo one can express the fidelity J- in terms of 



the surface free energy fs as = e ^ and therefore 
Xn{gi,gf) = dg^fs{gi,gf)- This expression clearly reveals 
that xi is the so-called excess internal energy while the 
fidelity susceptibility X2 is the excess specific heat [l^ of 
the d + 1-dimensional classical confined system. For the 
latter, it is known that X2 — \g — 5c I""" close to gc, 
where as is the excess specific heat exponent given by 
as = a -\- V where a and v are the corresponding 
standard bulk critical exponents of the classical c? -f 1- 
dimensional system. They satisfy the hyperscaling rela- 
tion a + i/{d+l) = 2 which yields immediately the scaling 
in Eq. ^ . In passing, we mention that mixed derivatives 
of In with respect to gf and gi are also expected to be 
characterized by critical exponents which generally differ 
from as even though, in most of the cases, they can be 



1^ 



expressed in terms of known bulk exponents 

Setting forth an argument based on the physical anal- 
ogy between critical phenomena in confined systems and 
quantum quenches, one can now show that P[W > Ai^o) 
close to the threshold Ai?o and to criticality 17 ~ is 
universal and determined by the asymptotic expression 
of the critical Casimir free energy fc- Close to a critical 
point, fc takes the scaling form 



fc{R) = R-''F{R/^+), 



(5) 



where ^+ coincides with the exponential correlation 
length of the system in the bulk disordered phase and 
F{x) is a universal scaling function which depends only 
on the surface and bulk universality classes and not on 
the detailed values of gi, gj, etc. While the large- i? 
decay of fc is oc R~'^ at criticality, away from the crit- 
ical point it is dictated by the asymptotic expansion of 
F{x^l). Generically, it takes the form 



F{x) = A{x) e-«^ + 



(6) 



where the power law A{x), q € M.'^, as well as F, de- 
pend, inter alia, on being above or below the bulk crit- 
ical point. For the quantum (classical) Ising model in 
d^ 1 (2) q e {1,2} where the correlation length 
= a\g — gd '^ is identified with the inverse of the mass 
TO of the lightest particle of the quantum model in the 
paramagnetic phase. Performing the analytic continua- 
tion Ri-^ it we obtain for G{t): 



Git) = e 



_ -iAEot-Lx2f^ 



(1 + A{tmt) e-"""* + ...), (7) 



which immediately shows that P{W) consists of a delta- 
function peak aA, w = W — AEq = and of a continuum 
for w > qm, where P{W) close to the threshold is deter- 
mined by A{x). It is important to stress at this point 
that while one would be naively tempted to associate the 
small-w behavior of P{W) to the large-t expansion of 
the generating function G{t) - which does not necessar- 
ily coincide with the one on the imaginary axis Eq. ([6]) - 
the correct procedure, instead, is the one described above 
and discussed in more detail on the example given below. 
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It is now useful to illustrate the ideas presented above 
in the specific case of the quantum Ising chain subject to 
a quench of the transverse field [l^ . The Hamiltonian is 



(8) 



where a^'^ are the spin operators at lattice site i (with 
lattice spacing a = 1), J (= 1 in what follows) is an 
overall energy scale, and g is the strength of the trans- 
verse field. The one-dimensional quantum Ising model 
is the prototypical, exactly solvable example of a quan- 
tum phase transition [isj . with a quantum critical point 
at gc — 1 separating a quantum paramagnetic phase at 
g > gc from a ferromagnetic one at g < gc- H{g) can 
be e asily analyzed after a Jordan- Wigncr transforma- 
tion |15| in terms of fermionic operators Ck of momen- 
tum k: H = 2J2k>a H Hk ^k, where = (cfc,cLj^, 
and Hk = [g — cosk)&z ~ (sinfc)(Tj,. This 2x2 matrix 
can be diagonalized with a rotation U{9) = cxp(— i^tTj;) 
of an angle 9k = l/2arctan[(sinfc)/((7 — cosfc)] which 
gives W{dk)HkU{dk) = sign{g-cosk)Ekaz, where Ek = 
[(g — cosfc)^ + sin^ fc]^/^. Accordingly, H takes the form 



H = J2 Ek hhk - 7-fc7lfc), 



(9) 



fe>0 



where we have conveniently introduced the quasi- 
particlcs jk^s in terms of the c/c's: 



7fc = (7^=,7lfc)^ = c''>'UH9k + Sk)-^k, 



(10) 



with Sk = for g — cos fc > and Sk = tt/2 otherwise. 

The statistics P{W) of the work W done during a 
quench of the transverse field g from gi to gj can be cal- 
culated by expressing |^'o(gi)) in terms of the eigenstates 
l*n(5/)> of H{gf). This can be done by using Eq. (HOl) 



which yields the relation jkigj) = e'^"''-" U^Aak) ^k{gi) 
between the operators jkigi) and ^k{gf) which allow 
the diagonalization of H{gi) and H(gf), respectively, 
where ASk = 5k[gf) - 5k{gi) and Aa^ = Ok{gf) - 
9k{gt) + A4. Since by definition 7fc(5i)l*o(gj)) = 
e-'''^*'= [cos(Aafe)7fc(5/)-isin(Aafe)7lj(7/)] \^o{gi)) - 
0, one can finally express |\I'o(.gi)) as 

|*o(g«)) =A/'e^'=>«**^"('^"'=)'''^(3-f)^-'=(f/)|^ro(5^))(ll) 



4fc>0 



In (l + tan'^ Aofe)] is a nor- 



where N = cxp[— | 

malization constant [16[ and the sum runs over the mo- 
menta in [0, tt] allowed by the finite extent of the system. 
Inserting now Eq. (dl]) into Eq. (P) one finds 



InG'(t) = -iAEot + ^ In 



A:>0 



1 + (tan Aafc)2e 



~2iEk{gf)t 



1 -I- tan"' Aak 



(12) 

According to Eq. ((TT|) the fidelity is given by 
-^(54,3/) = 1-^1 = exp[-Lfs{g^,gf)], where fs{gi,gf) = 



{2L)~^ Sfe>o [f + ^^^^ Aak] is nothing but the sur- 
face free energy of the classical two-dimensional (semi- 
infinite) Ising model. The successive derivatives with 
respect to gf of this expression give the first two sus- 
ceptibilities, analogous to the excess entropy and specific 
heat, e.g.. 



Xi 



^9ff = j^^M^(^k)dgfek{gf). 



(13) 



fc>0 



Close to criticality these sums are dominated by low- 
momentum modes and therefore they can be evaluated 
in the scaling limit from the outset [17[, which yields 



9k{g ^ 1) ^ ±(1/2) arctan[fc/TO(,9)] 



(14) 



for 9k at small k, where m(g) = |.g — 1|. Here we focus on 
critical quenches with ruf = rn{gf) <^ rn{gi) = rrii. For 
gi.f > 1, in Eq. jlS]) tan(Aafc) — fc/(2m/) for k < m/, 
~ 1 for rrif <^ k <^ rrii, whereas ~ mi /{2k) for <^ k. 
Taking into account Eq. (fT4)) one therefore finds 



Xi(.g/, .gi) = —l/(47r) ln(mi/m/) -I- regular terms. (15) 

(If the lattice spacing a, presently set to 1, is reinstated 
into the calculation, the prefactor of the logarithm be- 
comes (X a and is therefore not universal.) The same 
result is obtained for gi > 1 > gj , whereas xi has the 
opposite sign compared to Eq. ([15]) whenever gi < \. A 
further derivative with respect to gj gives the singular 
behavior of X2 — ^/{gf ^ 1) in agreement with the scal- 
ing law Eq. ([4]) in = 1 with the proper value = 1 in 
two dimensions. 

The final goal of our analysis is to calculate the statis- 
tics of the work W close to the threshold AEq and to 
criticality gi f :^ g^ focussing on the t-dependent part 
fc{gf,g^,t) = -l/LEfc>oln [l + (tanAafc)2e-2^^^(9/)*] 
of the exponent in the characteristic function. In order to 
extract its low-w asymptotics we consider as before the 
critical Casimir free energy fc{gf,gi,—iR), we extract 
its asymptotic expansion for large R and finally we con- 
tinue analytically an expansion of the form of Eq. ([7]). 
Since this asymptotics is again determined by the low- 
momentum modes it is convenient to take the scaling 
limit, which gives 

+ 00 



l + (tanAafc)2e-V'='+™/« , (16) 



fc = -l -In 
■' ' 2n 



where the 9k^s in Aak are given by Eq. (fT4| . Note that 
since our focus is on universal quantities, the expression 
above can be regularized at large momenta via a con- 
venient cutoff scheme. In particular, we imply that in 
doing the final analytic continuation R ^ U + tq we keep 
a small real part, which effectively cuts off the contribu- 
tion of ultraviolet modes in Eq. (fT6|). 

In order to classify the various possible \ow-w be- 
havior of P(VF), we consider the small-fc expansion of 
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ta,n(Aak) in Eq. (|16p . Below we summarize the main 
results of this analysis, the details of which will be 
presented elsewhere [18[. When the initial and final 
states are both either ferromagnetic or paramagnetic, 
the problem can be simplified expanding in the differ- 
ence gf — Qi, and tan(Aafe) ~ zLp-k/{2'mf) for gij ^ 1, 
where p± = [(m^ ± mf)/mi]^^. Calculating the asymp- 
totics of Eq. (|16p for x = 77i/i?^ 1 one obtains the scaling 
function in Eq. (O 



F{x > 1) ~ 



32w¥x 



(17) 



For a two-dimensional classical Ising model in a strip a 
critical Casimir free energy with the asymptotic depen- 
dence on a; = R/^+ ^ 1 given by Eq. (fT7|) with mf <^mi 
corresponds to the case of fully magnetized boundaries 
(effectively realizing the so-called (-I-, -I-) boundary condi- 
tions [l2j) and T <Tc, where T is the temperature and Tc 
its critical value in the bulk. By duality, this corresponds 
also to the case of effectively paramagnetic boundaries 
(realizing the so-called (O, O) boundary conditions [l2|) 
and T > Tc [l^ [13 ■ both these instances the bound- 
aries have the same paramagnetic/ferromagnetic charac- 
ter as the bulk. As expected, instead, F vanishes for 
rrii = TTif, as no quench has occurred. After the ana- 
lytic continuation, the characteristic function can be ex- 
panded as in Eq. ([7]) and one obtains for P{W) a peak 
at the threshold AEq plus a continuum starting at 2m, 
which describes pairs of quasi-particles. The correspond- 
ing edge singularity is universal and fully determined by 
the exponent of A{R), i.e., 



P{w)^S{w) + ^ ^ l^plJ i + (18 

4 V m/ 

where S = At: / L is the two-particle level spacing, w = 
W ~ /S.Eq, and § is the unit step function. Perturbation 
theory breaks down as soon as cither the final state is 
critical or the quench takes the system across the phase 
transition. In the first case, if = 1 (or gi ~ 1), 
tan(Aa/j) ~ 1, while 



tan(AQ!fc) ~ ±2p_|-TOyfc 



(19) 



for gi^l^gf In both cases the behavior of P at small 
w is significantly different from the one obtained above. 
Indeed, while for quenches at criticality one obtains a 
scaling function F{x = 0) = — 7r/48, consistent with the 
results of conformal field theory [l3| , for quenches across 
the quantum critical point and a; 3> 1 one has 



Fix) 



-p+xc 



2 S. - 

zp_^_x'e 



Eq. ([20|) with mj <^ is observed in the critical 
Casimir effect for T > Tc and (-I-, -|-) boundary con- 
ditions or, by duality, for T < Tc and {O, O) bound- 
ary conditions, i.e., whenever the boundaries and the 
bulk have opposite paramagnetic/ferromagnetic charac- 
ters. Comparing the two cases corresponding to Eqs. pT)) 
and ((20| for nif <^ m^, we note that the associated re- 
sults for F{x) are consistent with a change in the effec- 
tive boundary conditions upon changing the ferro mag - 
netic/paramagnetic character of the initial state [l8, 21 



Taking now the analytic continuation of the previous re- 
sults we obtain 



P{w) (X S{w) 



125 



(21) 



for quenches to the quantum critical point, i.e., the con- 
tinuum starts exactly at the edge, as expected from 
having gapless excitations in the final state, while for 
quenches across the quantum critical point we have 

2n 

P{w) (X S{w) + —(2mfp^)5{w — ruf) 



'd{w—2m) 



(2p+)' 



w — 2raf 



(22) 



The emergence of the delta- function peaks at mf,2mf, 
etc.. is a non-perturbative phenomenon consequence of 
the kinematic pole in Eq. ([T9l) . Notice that while the 
edge singularity at 2m appears to be not integrable, the 
presence of a cutoff at w — 2m / « 1/L is implied by the 
finite size of the system. 

Summarizing, we established a connection between the 
statistics of the work done during a quantum quench and 
the critical Casimir effect. This relation allowed us to 
demonstrate how universal features emerge in this non- 
equilibrium evolution depending on the initial conditions. 
We expect that the extension of this analysis to differ- 
ent systems and protocols will help elucidating the role 
of universality in the dynamics of closed and strongly 
correlated quantum systems. 
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(20) 



where we kept also the leading term in exp(— 2a;), which 
controls the edge singularity associated to the two- 
particle continuum. The behavior corresponding to 
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